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We investigate theoretical decoherence effects of the motional degrees of freedom of a single 
trapped atomic/ionic electronically coded qubit. For single bit rotations from a resonant running 
^~>) , wave laser field excitation, we found the achievable fidelity to be determined by a single parameter 

' characterized by the motional states. Our quantitative results provide a useful realistic view for 

\ current experimental efforts in quantum information and computing. 
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Since the pioneering work by Shor on efficient prime 
factorization with a quantum computer in 1994, we have 
witnessed an explosive growth of interests in quantum 
information and computing. Although still largely a the- 
^ ' oretical field, solid progress in experimental efforts have 
. been made within the last few years. Most notably, pure 
■ state based quantum gate implementation ||^, demon- 
, strations of quantum teleportation Q, and GHZ state 
L.J ' synthesis have stimulated more vigorous experimen- 
' tal efforts. 

' Two of the most interesting proposals for potentially 
■""r!; ' large scale quantum computing were suggested by the 
same Innsbruck group, based on trapped ions |^ and 
cavity QED with atoms Various implementations 
of these ideas are actively pursued in many experiments 
around the world. In their original analysis as presented 
^ in dJ^, individual qubits are coded in electronic de- 
O^' grees states of atoms/ions, and coherent evolution of the 
^ ] system state requires the qubits to be in selected mo- 
• . tional pure states. Experimentally, one needs to attain 
/\ ' the strong binding limit j7| and cooling to the motional 
^ , ground state is also required As is well known, these 
- - ' limits are difficult to maintain due to various decoher- 
ence processes [|| which heat up the motional degrees 
of freedom. Further more, maintaining motional ground 
state becomes problematic when strong confinement is 
not satisfied. Several ideas were proposed recently for 
computing with 'hot' qubits | |To[ |. 

This paper attempts to provide quantitative answers 
to motional effects (ME) on electronically encoded quan- 
tum states It is the first step towards a thorough 
investigation of the ME. The paper is organized as fol- 
lows. First our model is presented. We then discuss 
analytically the decoherence of the ME for an unknown 
electronic encoded qubit. Finally we present numerical 
results to support our understanding. In forthcoming pa- 
pers, we will study the decoherence due to ME on multi- 
qubit entanglement creation, e.g, the effect on the condi- 
tional logic operation CNOT. 

We consider a single harmonically bound two state 
atom described by the Hamiltonian [p^-p^ 

H= V fi{n^,ujl + UyUjl + n,uf)\g, n){g, n\ 



rh 

+ inf]ie''^^*^77^;„-,(fci)|5,n)(e,m| + /i.c. . (1) 

\g,n) = \g)\n)g (|e,m) = |e)|m)e) denotes number state 
in the ground (excited) trap with frequencies are tff^^ y ^ 
(wf). ojeg is the electronic transition frequency. 17l is 
the Rabi frequency of the plane wave laser field. The 
motional dipole moments are the familiar Franck-Condon 
factor (l| 

VnrfXkL) ^ {9,n\e-'^--%,m). (2) 

Radiative coupling to the vacuum reservoir of the atom 
will not be included here as the effect of the resulting 
spontaneous emission on the qubit decoherence has been 
studied and is well understood [5|J^, p^ , |l5| . Our model 
can also be viewed as between two ground states in a 
three level A-type off resonant Raman system. In such a 
case, VLl ^ ripWg/SL are the two photon effective Rabi 
frequency and fc/, = kp — ks- The indices P and S denote 
the dipole connected pump and stokes transitions, and 
(5l is the (large) detuning from the eliminated far off- 
resonant excited state. 

Physical models similar to Eq. (0) have been stud- 
ied under different context before Rfzi ^^lJl^Jl^ , usu- 
ally, within the Lamb-Dicke limit (LDL) when ME be- 
come considerably simplified. In terms of the trap width 
af^ = ^/hj2MLof^ , the LDL corresponds to kLaf <C 1. 
This requires the atom/ion to be confined less than the 
wavelength a^'"^ ^ Al and is equivalent to require hcuf to 
be much larger than the recoil energy Ep — h^k\/2M <C 
fiLo^'^ ■ For an effective two state system reduced from 
a near resonant three level A-type configuration, LDL 
is easily satisfied with co-propagating pump and Stokes 
fields when kp ~ ks- 

In this study we investigate ME for general cases not 
in the LDL. Such studies will provide much needed theo- 
retical clarification as trapped atoms/ions are among the 
"hottest" qubit candidates in many experimental efforts. 
We note Franck-Condon factors (^ satisfy 

^ ^ [^/nm(fcL)] Vitrn' {kL) — ^mm' y 
n 

[Vnin{kL)\*Vn'm{kL) = Sflft' , (3) 

which allows the introduction of a complete and or- 
thonormal basis 

\^)p^^V*vn'{^L)\n)g. (4) 
n' 

Physically it corresponds to the motional wave packet 
of a photon absorption from ground state \n)g. Mathe- 
matically, it is the number coherent state basis |a,n) = 



D{a)\n), with the displacement operator D{a) = 
^ab+-a'b_ ^Q^g ^. ^ _^ ^^^Yi bl (6f+) the 

annihilation (creation) operator for|ni)g, therefore Eq. 

can be rewritten as \nx)p — \ikLa^,nx)g, exactly rep- 
resenting wave-packets corresponding to excitation from 
different ground trapping states. As will become clear 
later, the coherent Rabi coupling between the ground 
and excited state manifolds can also be decomposed into 
paired sets {\g)\n)g, \e)\ri)p}. 
With the inverse relation 

\m)e = V'??mn(fcL)|"^)p, (5) 



we can transform Eq. into the \n)p basis. Denote 

\ep,n) = |e)|n)p, we obtain 

?iWeg I e, m)(e , m I = ^ hujeg \ep,n){ep,n\, 

rh fi 

^miTu^l\m)ee{m\ = J^i^SrV + E2n')\n)pp{n\. (6) 

mji fin' 

We found SS,, terms couple nearest neighbors, i.e. states 
with Hi — n'^ ± 1 [rij^i = n'-), while E9^, terms couple 
states with fi — fi' and Ui = n[±2 (nj^i = n'j). E9^, 
terms become diagonal along i— axis whenever dojf 7^ 0. 

We will focus on the case u)f = uof {i — x^y,z) in the 
present paper. This is typical for the ion trap system and 
can also be arranged for optical dipole traps . For the 
plane wave excitation along the x-axis, the ME along the 
y and z directions are unperturbed. Our Hamiltonian 
Eq. (|l|) can simplifies to a one dimensional model 

H= ^n^huj^\g,n^){g,n^\ 

+ '^h{nxuj^ - AL)\ep,nx){ep,nx\ 
+i{kLax) ^ hcu^y/n-j; + l|ep, + l)(ep, n^l + h.c. 
+ '^^^L^\g,nx){ep,nx\ + h.c, (7) 
where transformation to the interaction picture by 



U{t) = exp -inLjLt^\ep,n^){ep,n^\ (8) 



has also been made. The detuning is HAl ~ hojL — hLj^g — 

^jj- , including the recoil shift. 

This Hamiltonian can be graphically illustrated as in 
Figure |l|. The paired ladder structure resembles the fa- 
miliar motional state ladders in an ion trap [Q. However, 
in Fig. |l|, the excited states are the wave-packet basis 
states Eq. (^). The nearest neighbor coupling is not due 
to the LDL approximation. Denote H = Hq+Hi with Hi 



the nearest neighbor couphng term in the excited state 
[the third line of Eq. (|^)], Hq becomes 



Hn = h 



fir 

-f- "X^; 



2 

3 - A, 



(9) 



where each 2x2 block describes the Rabi oscillation be- 
tween paired states {\g)\nx)g,\e)\nx)p} with exactly the 
same Rabi frequency = i/ri^ + A|. There are no dif- 
ferential detunings between different pairs either. With 
wave function coefficients {C,^^ , C^^}, the 2x2 oscillation 
is described by 

Cl(r)=e-»^-^- [C1(0) cos 0-zQ^(O) sine] , (10) 

when Ai = 0. The pulse area is 9{t) = ^ JJ" V,{t)dt. Hq 
describes the coherent evolution between paired states 
with a time scale given ft. Hi, on the other hand couples 
nearest neighbors of excited motional wave pack states 
and can causes decoherence of an electronically coded 
qubit. Its time scale is determined by several factors in- 
cluding the trap frequency wf, Lamb Dicke parameter 
k^a^, and the highest motional state number n^'^^. As- 
suming an electronically coded unknown qubit 

|V'(0))=a|5)+/?|e), (11) 

(normalization jap -I- |/3p = 1), an arbitrary single bit 
rotation is achieved through a multiplication of p9| ] 



s'^i \ / cos 61 sin \ / e'l^ 
e-^^i j \ -sine cose j \ 



^1 ll^ 72 7 and (j) are parameters. With resonant Rabi 
coupling, our solution Eq. ( |l0|) achieves the important 6 
rotation corresponds to e^*"*"^ = — i, and e^^^^ — i. Ideally 
one hopes to arrive at the target state 

|'(/'(r))T= (acos^ — ij3sm9)\g) + {(3cos9 — iasin6')|e), 

with the density matrix, 

PT(r)= Ij\g){g\ + (1 - /J)|e)(e| + (/Jjg)(e| + h.c), (12) 
where 

/J= lap cos^ 9+\(3\^ sin^ 9 + i{a(3* - c.c.) sin 9 cos 6*, 
/J^= a(3* cos^9 + Pa* sin^ 6* -f i(|ap - sin9 cos9. 

Due to ME the electronic qubit Eq. ( [ll| ) does not 
remain decoupled from the motional degrees of freedom. 
In general, it evolves within the much larger Hilbert space 
containing motional states. We now study two concrete 
examples of decoherence assuming the initial qubit in 



the enlarged Hilbert space reproduces the density matrix 



Eq. (12), i.e. resembles a perfect qubit to the innocent 
bystanders unaware of the motional degrees of freedom. 
First, we consider 

|i^(0))tot= («|<?) + P\e)) ® l^(0))c„., (13) 

with an initial pure motional state |'!/'(0))cm — 
Sn, c„^|na;)g {Y.n^ |c«xP = !)■ Upon tracing the mo- 
tional degrees of freedom, p(0) — |i/'(0))(V'(0)| is correctly 
reproduced. By rewriting |V'(0))tot as 

Vn^n'^Cn'^ P\S) \ nx) p ) , 



we can analytically evolve this state with Hq to obtain 
p(r)= Ig\g){g\ + (1 - /,)|e)(e| + (/,,|g)(e| + h.c), (14) 



with 



Ig= lap cos^ sin^ 6 + i(a(3*ri* - c.c.) cos 6 sin 

/,e=*(|ap-|/3nisin20 ^ e-^^"^-''^^"^^c„,c;^,, 

+a/3*cos2^ ^ e-'("-^^)-^-c„^^77;,,c*,77,^„. 



+a*/3sin2 6' ^ e-'("==-«-)"'^ ^ 



,gx 



where we have defined the parameter 



We see the evolution by Eq. (^^ will in general not 
reproduce the intended density matrix Eq. (03) because 
of rj. Hi term is the other reason for incomplete control 
although its effects (when fl^ 3> lo^) can be minimized 
by employing a fast pulse with w|r <^ 1. Within such a 
limit, or when ajfr — 2-k, we obtain 

/ge= i{\a\^ - |/3|^)sin6'cos6l?7(fcL) 

+a/3*cos2 61 + a*/3sin2 61 ?7(2A:i). (16) 

The necessarily condition for attending a perfect fidelity 
of the single bit rotation is then rjikL) = 1, which can 
be approximately satisfied in the LDL when fciag <C or 
in the A-type Raman systems with co-propagating pump 
and Stokes fields. As a second example, we consider the 
case of a thermal motional state 

Ptot(0)= |7/;(0))(^(0)|®pcm(0), 

J^x 

This is the limiting case of an ensemble average of Eq. 
( |l3| ) with c„^ = ^yO™ e~*'^"=i= and 4>n^ a uniform ran- 
dom number e [0, 27r). The dynamics due to Hq can 



be evolved analytically and the same density matrix Eq. 
( p^ ) is obtained. After averaging over {</>„}, we obtain 



exp 



1,, ,o , fl huis 



(18) 



In the low temperature limit when k^T < tiuj^, rj be- 
comes 1 as long as LDL /clo^ <C 1 is satisfied. At high 
temperatures when ksT » {kLO,^)'^ <^ hcj^/kBT 

needs to satisfied for ij close to 1. 

We now discuss the numerical solutions. Expand the 
total wave-function as 

|^(i)>tot = ^[cL(i)l.9)l«.)ff + <Jt)|e)|n,)p], (19) 



we have solved the Schrodinger equation including both 
Ho and Hi. The transformation Eq. (Q) greatly reduces 
the motional Hilbert space dimension. The perceived 



fidelity for the electronic coded qubit Eq. (11) under 
transformation Eq. ( p^ ) is 

T= Ty[pt{t)p{t)]. (20) 

We take a = (3 — 1/V2 as an example to illustrate 
our numerical results since similar/better fidelities are 
obtained with other choices. In Figure ^ we compare fi- 
delities under arbitrary 9{t) rotations for two different 
pure states. Acceptable fidelities are obtained only for 
fciof < 0.3 . In general larger Ql/^^x ratios also im- 
proved fidelity although it saturates around VIl/oj^ ~ 
100. Noticeable improvements are also recorded for nar- 
rower distributions in |c„p, e.g. in Fig. |^ initial state 
with c„^ ~ (5„^,o produced better fidelity. This is a di- 
rect reflection of dephasing among different motional pair 
states because of their different time scales from Hq and 
Hi . The oscillatory behavior is due to dephasing caused 
by the Rabi oscillation between motional paired states. 
For comparison, we note that meaningful single bit ro- 
tations need to achieve a fidelity of 1/2, the lower limit 
from a random (uncontrolled) sampling of final states. 

Finally we compare with thermal states for several dif- 
ferent values of ksT/huj^.. Surprisingly, we found the 
fidelity for an initial motional thermal state is always 
higher than its corresponding pure state. In the tem- 
perature regime considered we found acceptable results 
as long as LDL is maintained. In Fig. we have used 
kLa3 = 0.1 and = 100 (wf). 

In conclusion, we have performed detailed theoreti- 
cal studies of the decoherence of an electronically coded 
atom/ion qubit due to ME. By introducing a wave packet 
basis in the excited state, we were able to perform con- 
siderably cleaner analysis to simplify the ME. We found 
that a single parameter rj measures the achievable fidelity 
of arbitrary single bit rotations. We performed numerical 
calculations which demonstrates our understanding and 



provided quantitative limits for experiments: the LDL 
is always required to maintain a high fidelity for arbi- 
trary single bit rotations. We also found that a pure mo- 
tional state is not necessarily preferred although a qubit 
with an initial ground motional state does give rise to the 
highest recorded fidelity. In actual experimental imple- 
mentations, a large 17^ is also needed to assure negligible 
motional dephasing during r. One can always wait for a 
period Stt/w^ for subsequent single bit operations since 
the motional wave function then always rephases to its 
initial state. This study will also shed light on devis- 
ing schemes for overcoming ME decoherence and error 
corrections in trapped atomic/ionic qubits. 
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FIG. 1. The paired ladders of a trapped two state 
atom/ion. Solid curve arrow heads denote Rabi oscillations 
between paired states, while the dotted curve arrow heads 
denote nearest neighbor motional coupling. 
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FIG. 2. The fidelity for single bit rotations with fc^af — 0.1 
(dashed lines), 0.3 (solid lines), and 1.0 (dot dashed lines). 
The higher and lower fidelity sets are for initial motional state 
Cn^ = Sn^o and c„^ = (2(5u^o + V^Sn^i + 5n^2)/V7 respec- 
tively. Ql = WO {u^). 
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FIG. 3. The fidelity for single bit rotations with 
kBT/hui^ = 1 (dashed lines), 3 (solid lines), and 10 (dot 
dashed lines). The higher and lower fidelity sets are for ini- 
tial motioal thermal state p™ and pure state c„^ — \/p^ 
respectively. 



